Detectable inertial effects on Brownian transport through narrow pores 
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We investigate the transport of suspended Brownian particles dc driven along corrugated narrow 
channels in a regime of finite damping. We demonstrate that inertial corrections cannot be neglected 
as long as the width of the channel bottlenecks is smaller than an appropriate particle diffusion 
length, which depends on both, the temperature and the strength of the dc drive. Therefore, 
transport through sufficiently narrow constrictions turns out to be sensitive to the viscosity of the 
suspension fluid. Applications to colloidal systems are discussed. 
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Brownian transport in narrow corrugated channels is 
a topic of potential applications to both natural [H [2] 
and artificial devices . Depending on the amplitude 
and geometry of the wall modulation, corrugated chan- 
nels fall within two distinct categories: (i) smoothly 
corrugated channels, typically modeled as quasi one- 
dimensional (ID) periodic channels with axial symmetry 
and unit cells delimited by smooth, narrow bottlenecks, 
also called pores [iUl2|: (ii) compartmentalized channels 
[I3HT6| . formed by identical compartments separated by 
thin dividing walls and connected by narrow pores cen- 
tered around the channel axis. Brownian transport in 
such sharply-corrugated channels must be treated as an 
irreducible two or three dimensional diffusion problem 
[T7] . More importantly, for both categories of corrugated 
channels most analytical results only apply under the 
condition of very narrow pores [2 [3 [17] . 

Corrugated channels are often used to model transport 
of dilute mixtures of small particles (e.g., biomolecules, 
colloids or magnetic vortices) in confined geometries 
Each particle is subjected to thermal fluctuations with 
temperature T and large viscous damping 7, and a ho- 
mogeneous constant bias of strength F parallel to the 
channel axis. Such a dc drive is applied by coupling the 
particle to an external field (e.g., by attaching a dielectric 
or magnetic dipole, or a magnetic flux to the particle), 
without inducing drag effects on the suspension fluid. In- 
terparticle and hydrodynamic interactions are thus ig- 
nored for simplicity (these assumption are discussed in 
Ref. [4). 

In this paper we investigate the relevance of inertial 
effects due to the viscosity of the suspended particle. 
As is often the case with biological (and most artiflcial) 
suspensions the Brownian particle dynamics in the 
bulk can be regarded as overdamped. This corresponds 
to (i) formally setting the mass of the particle to zero, 
TO = 0, or, equivalently, to make the friction strength 
7 tend to infinity, and (ii) assuming F smaller than the 
thermal force Fq = ^y/kT/m (Smoluchowski approxima- 
tion) [18, . The current literature on corrugated channels 
invariably assumes such an overdamped limit. But how 
large is an "infinite" 7 (or how small can be a "zero" to)? 
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FIG. 1: (Color online) Rescaled mobility 7^ versus the 
rescaled 7 unit, rj/jr, of Eq. (|5| in the corrugated channel 
of Eq. (pl for small drive, F/ Ft ~ 0.1, and for different pore 
sizes, A/j/L, and aspect ratios, r = xl/vl- The relevant scal- 
ing parameters are Ft = kT/A and 7t = V mkT/A, Eq. (jsj. 
The dotted line is the fitting law of Eq. ([5|. The actual pa- 
rameter values used in our simulations are m = kT = yL ~ 1. 
Inset: corrugated channel with w{x) as in Eq. ([2|. 



The answer, of course, depends on the geometry of the 
channel. Our main conclusion is that the overdamped dy- 
namics assumption for Brownian diffusion through pores 
of width A sub jected to a homogeneous drive F, applies 
only for 7 3> V mkT/A, and 7 ^ y/mF/A, irrespective 
of the degree of corrugation. This means that the inertial 
effects cannot be neglected as long as the Brow nian d if- 
fusion is spatially correlated on a length (It = V mkT/j 
at small dc drive, or = mFjj^ at large dc drive) of 
the order of, or larger than, the pore width A. Therefore, 
for sufficiently narrow pores or sufficiently large drives, 
inertia always comes into play by enhancing the blocking 
action of the channel bottlenecks. Thus, the condition of 
vanishingly narrow pores, A 0, assumed in most ana- 
lytical studies, can be inconsistent with the assumption 
of overdamped diffusion in the spirit of the Smoluchowski 
approximation |18) . Finally, we also discuss applications 
to transport in colloidal systems. 

Let us consider a point-like Brownian particle diffusing 
in a two dimensional (2D) suspension fluid contained in 
a periodic channel with unit cell cci x y^, as illustrated in 
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FIG. 2: (Color online) Rescaled mobility, 7/1, in a corrugated 
channel with r = 1, A/yL ~ 0.1, and (a) vs. the damping 
7/7T for diflerent F; (b) vs. F/Ft for various 7. 7t, Ft 
and the actual simulation parameters are as in Fig. [l] In- 
set (a): data from the main panel vs. the rescaled damping 
constant, 7/7F, with 7f — \/mF/A, Eq. the raising 

branch of the universal curve for large drives is fitted by the 
power law (7/7f)°' with a = 1.4. Inset (b): 7/i vs. A/yL for 
7/7T = 0.4 and two values of F. The dashed lines in (b) rep- 
resent the large-drive power laws {F/Ft)~'^^^ (main panel) 
and (A/yi)"/^ (inset), both with a = 1.4; for F/Ft = 1, M 
is proportional to A (inset). 



Fig. [I] (inset). The particle is subjected to a homogeneous 
force F = Fcx oriented along the x-axis. The damped 
Brownian particle obeys the 2D Langevin equation. 



■m(fr/df = -7 dr/dt + F + y^2^ ^{t), (1) 

where f— {x,y). The random forces ^(t) = {S,x{t) , £,y{t)) 
are zero-mean, white Gaussian noises with autocorrela- 
tion functions {^i{t)£,j{t')) = SijS{t — t'), with i,j — x,y. 
The symmetric walls of the corrugated channel have been 
modeled by the sinusoidal function ±'w{x) (Fig.jlj inset), 

w{x) = (1/4) [(2/i + A) - {yL - A) cos(27ra;/a;i)]. (2) 

With this choice for w{x) we made explicit contact with 
the current literature on entropic channels [4Hl0]. How- 
ever, similar results were obtained for sharper corruga- 
tion profiles, e.g., for compartmentalized channels |13H16| 
(not reported here). We numerically integrated Eq. ([!]) 
by a Milstein algorithm [19] with time-step, At, ranging 
from 10^^ down to 10~^, as F was increased. For each 
run, At was set small enough for the output to be in- 
dependent of it. The stochastic averages reported below 
were obtained as ensemble averages over 10^ trajectories 
with random initial conditions; transient effects were es- 
timated and subtracted. 



Two quantifiers have been used to characterize the 
Brownian transport in such sinusoidally corrugated chan- 
nel: 

(i) mobility. The response of a Brownian particle in 
a channel subjected to a dc drive, F, oriented along the 
X-axis, is expressed by the mobility, n{F) = {v{F))/F, 
where {v} = {x{F)) = limt^oo[(a;(t)) - x{0)]/t. The 
function ^i{F) increases from a relatively small value for 
F = 0, i.e. fi{F = 0) = up to the free-particle limit, 
7Aioo — 1, for F — >■ cxD [lOj. We recall that in the bulk, a 
free particle drifts with speed F/7. 

(ii) diffusivity. As a Brownian particle is driven across 
a periodic array of bottlenecks or compartment pores, its 
diffusivity, D{F) = \\mt^^[{x'^ {t)) - {x{t))'^]/2t, picks 
up a distinct F dependence. In corrugated channels with 
smooth bottlenecks, for F 00 the function D{F) ap- 
proaches the free or bulk diffusion limit, D{oo) ~ Dq = 
kT/j, after going through an excess diffusion peak cen- 
tered around an intermediate (temperature dependent 
[lOj ) value of the drive. Such a peak signals the depinning 
of the particle from the barrier array of w{x) associated 
with the channel bottlenecks [[20J. 

Recall that, in the absence of external drives and for 
any value of the damping constant, Einstein's relation, 
7/io — D{F — 0)/Dq, establishes the dependence of the 
transport parameters on the temperature and the chan- 
nel compartment geometry under equilibrium conditions. 

Inertial effects in corrugated channels become appar- 
ent both for small 7 and for large F. By inspecting Figs. 
[1] and [2] we immediately realize that (for small 7) inertia 
tends to suppress the particle mobility through the chan- 
nel bottlenecks. At large F, when plotted versus 7 [Fig. 
[2ja)], the rescaled mobility approaches unity for 7 — > 00, 
as expected in the Smoluchowski approximation [lOj , but 
drops to zero in the underdamped limit, 7 — ?> 0. More 
remarkably, the resulting 7// curves shift to higher 7 on 
increasing F (main panel). On expressing 7 in units of 
7t = \J mkT I A, see Eq. ([3]) below, the curves for large 
drives tend to collapse on a universal curve well fitted 
by the power law (7/7T)" with a = 1.4 (inset). Corre- 
spondingly, in Fig. ^b) the mobility grows hke A"/2 at 
large F (inset) and decays like F~"/^ for small A (main 
panel). 

The power law, 7/.* cx (7/7T)", introduced here is only 
a convenient fit of the rescaled mobility function, even if 
it holds for two or more decades of 7/7T. The analytical 
form of that function remains to be determined. Devia- 
tions from the fitted power law, /i oc A"/^, for the mo- 
bihty at very small A [Fig. [2jb), inset], suggest that the 
fitting exponent, a, slightly depends on A, with a — >■ 2 
in the limit A — (not shown). [Note that, on the con- 
trary, the power law 7/i oc F~"/^ works throughout the 
entire F range explored in Fig. [2jb).] 

The dependence of the rescaled mobility on the sys- 
tem parameters in the underdamped limit is further il- 
lustrated in Fig. [1] where at low 7 and for vanishingly 
small drives, the mobility grows proportional to the as- 
pect ratio r = xl/j/l of the channel unit cell and the 
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pore cross section A, i.e., scales with the dimensionless 
quantity r^/^x [see also the data set for F/Ft = 1 in 
the inset of Fig. [gb)]. 

Deviations of the diffusivity data from the expected 
overdamped behavior are even more prominent. As 
shown in Fig. [sj at large 7 the curves D{F) approach 
the horizontal asymptote D{F) — Dq, as expected pUj . 
However, beyond a certain value of F, seemingly pro- 
portional to 7^ (see figure inset), these curves abruptly 
depart from their horizontal asymptote. In the under- 
damped limit, the F dependence of the diffusivity bears 
no resemblance with the typical overdamped behavior. 
At low 7, all D{F) data sets collapse to a unique curve 
[Fig. [sj inset] , which tends to a value smaller than Dq for 
F — > 0, and diverges for F — ?> cxo, like F^ with /3 ~ 1. 
Such power law holds for large 7, as well, though in the 
large F domain, only. Indeed, for exceedingly large F, 
all D{F) curves seem to eventually approach a unique 
asymptote, irrespective of 7. 

By comparing the plots of Figs. [T][3] we conjecture 
that corrections due to inertia become significant in two 
regimes: 

(i) at low drives for 
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— V mkT / A 



(3) 



This characteristic damping was used to rescale the mo- 
bility data in Fig. [l] [see also Fig. ^h), inset]; moreover, 
in Fig. [3] for 7 < 7^ the diffusivity becomes a monotonic 
function of F with no plateau around Dq. The physical 
meaning of jt is simple. For 7 < 7^ the thermal length 
It ~ V mkT/ J grows larger than the width of the pores, 
A, so that the Brownian particle cannot reach the nor- 
mal diffusion regime, implicit in Einstein's relation, be- 
fore bouncing off the pore walls. As a consequence, the 
Smoluchowski approximation fails in the vicinity of the 
bottlenecks. 

Replacing 7 with 77^ in the bulk quantities Dq and 
Fq yields, respectively, Dt — kT/jT and Ft ~ kT/A. 
These are the 7-independent rescaling factors introduced 
in Figs. l][3 to characterize the inertia effects of the pore 
constrictions; 
(ii) at high drives, for 
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< 7f = y/niF/A 



(4) 



In the presence of a strong dc drive, the condition 7 ^ 7t 
does not suffices to ensure normal diffusion; the addi- 
tional condition A > Ip is required. Here, Ip — mF/j'^ 
represents the ballistic length of a driven-damped parti- 
cle, which is an estimate of the bouncing amplitude of 
a driven particle against the bottleneck. Upon increas- 
ing F at constant 7, Ip eventually grows larger than A 
and inertia comes into play. This mechanism is clearly 
responsible for the abrupt increasing branches of D{F) 
in Fig. d 

An analytical derivation of the transport quantifiers in 
the presence of strong inertial effects (low 7 and/or large 




FIG. 3: (Color online) Rescaled diffusivity, D/Do, vs. F/Fq 
(main panel) and D/Dt vs. F/Ft (inset) in the corrugated 
channel of Eq. ([2| with r = 1, A/y^ = 0.1, and different 
7. The actual simulation parameters are m = kT = yL = 1. 
The scaling parameters, Dt = kT/'yT, Fq = ^ \J kT / m, and 
Ft = kT/A, are introduced in the text. The solid line in the 
inset is the heuristic power law of Eq. Q . 



F) proved a difficult task. However, the following phe- 
nomenological arguments can explain well our numerical 
findings. 

Let us first consider the rescaled mobility at low drives. 
For F = the transport quantifiers 7/io and D{0) can 
be formally expressed in terms of the mean exit time, fg, 
of the Brownian particle out of a single compartment, 
namely, D{0) = a;|/4fe and ^iq = D{0)/kT (Einstein's 
relation). An analytical expression for fe as a function 
of the compartment geometry is only available in the 
Smoluchowski approximation |17j . In the low damping 
regime, we only have a rough estimate of fg. For F ^ 
the particle bounces off a compartment wall with rate 
(2/7r)(2-y/fcT/m/xi) (attack frequency), but only a frac- 
tion A/dt of such collisions leads to pore crossing. As a 
consequence, fe ~ A,XTyL\/ m/kT /ix A and 
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which closely reproduces the numerical data of Fig. [T] for 
F/Ft = 1 with mobility ^ oc A. 

The divergence of the diffusivity with the drive shown 
in Fig. [sj D{F) oc F^ with /3 = 1, is surprising since 
in the Smoluchowski approximation D{F) always tends 
to its bulk value, Dq. A heuristic explanation for this 
asymptotic power law runs as follows. The dispersion 
of a strongly driven particle at low damping is due to 
its bouncing back and forth inside the channel cells with 
speed ztv and v — F/j. Therefore, D{F) = w^fe/4, 
where fe is the particle mean exit time through one of 
the pores of a channel cell, fg = {7t/2){xl/v){'j/jt)- The 
time constant fe is estimated here as the cell crossing 
time, xl/v, multiplied by the geometric factor 7r/2 [due 
to the sinusoidal profile w{x), also used to derive Eq. 
([5])] and the success probability, 7/7T, for the particle to 
actually cross the bottleneck during a sequence of corre- 
lated bounces extending over the relaxation time m/7. 
One thus concludes that D{F) = {-k /%){xlAF) / \/mkT , 
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or, in rescaled units of 7t and I?t, 
D{F) _ ttxlF 



kT 



(6) 



in rather good agreement with our simulation data (see 
inset of Fig. pi). Lowering the temperature, for small 
damping diverges like T~^l'^^ which means that 

diffusion is dominated by a (chaotic) mechanism of bal- 
listic collisions. 

The main result of this work is that for real physical 
suspensions flowing through compartmentalized geome- 
tries, both in biological and artificial systems, the limit 
of vanishingly small pore size becomes extremely sensi- 
tive to the finite viscosity of the suspension fluid. This 
is not a purely academic conclusion. Inertial effects can 
be directly observed, for instance, in a dilute solution of 
colloidal particles driven across a porous membrane or an 
artificial sieve [5T]. To be specific, let us consider spher- 
ical polystyrene beads [22] of radius rg G [0.1, l]/im, sus- 
pended in a low viscous medium with, say, r\ — 8.9 x 10^* 
Pa-s (water), 6.04 x 10^* Pa-s (benzene) or 3.6 x 10^** Pa-s 
(acetone). For the typical mass density of polystyrene, 
200 kg/m'^, the mass of a spherical bead of radius O.l/im 
is m ~ 8.4 X 10^^^ kg. The corresponding rescaled 
damping constant 7/m, Eq. (IT]), can be determined by 
making use of Stokes' law, "ffm = 67Tr]ro/m, namely 
^/m= 1.9x 10"'^ ps~^ (water), 1.3x10""^ ps~^ (benzene) 
and 7.6 x 10"'* ps~* (acetone). 

At room temperature, T = 300 K, the thermal energy 
unit is kT = 4.14 x 10"^* kgm^/s^. If the pore radius 
is, say, 20% larger than rg, then the effective width A for 
the fine-radius particle to go through the pore, becomes 
A = A — 2ro = 0.02 /im and the rescaled zero-drive 

damping threshold jx/'m = \J kT /mh? is 3.6 x 10~^ 

ps~*. Moreover, a colloidal particle can carry an effec- 
tive charge of up to 10^|e|, e denoting the electron charge. 
Therefore, in the presence of an external electric potential 
of 10'^ — 10"* V/m, a colloidal particle experiences a pulling 
force of the order of 1—10 pN. The corresponding non- 



scaled large-drive damping threshold ^p/m — y F/mA 
then falls within the range (1 - 3)xl0~^ ps~*. In con- 



clusion, both 7t and are 10 to 100 times lower than 
the estimated particle damping constant 7/rn, so that 
the diffusion of such very small colloidal particle through 
a narrow pore (as for most biomolecules) can be safely 
described assuming an overdamped dynamics. 

However, by increasing the particle size, inertial effects 
come eventually into play. By the quantitative estimates 
reported above, for a colloidal particle 10 times larger in 
size, tq = 1/im, one obtains 7/771 — 1.9 x 10^^ ps^* (wa- 
ter), 1.3 X 10^^ ps^* (benzene) and 7.6 x 10^^ ps^* (ace- 
tone). This means that the damping constant 7 for such 
bigger colloidal particles falls between the corresponding 
thresholds 77- and 'jp, as now 7t/to = 2.2 x 10~^ ps~* 
and ^p/m. = (0.4 - 1) xlO~^ ps~*. Moreover, for sus- 
pended particles of even larger radius or made of mate- 
rials with higher mass density, the condition 7 < 7t,7_f 
is also experimentally achievable. As a result, for appro- 
priate values of the suspension fluid viscosity, a narrow 
corrugated channel can be used to efficiently segregate 
the constituents of a dispersed colloidal solution. 

In conclusion, the experimental demonstration of iner- 
tial effects on Brownian transport through narrow pores 
is easily accessible by manipulating artificial particles of 
moderate size by means of well established experimen- 
tal techniques. For too small particles, like biological 
molecules, detecting such effects might require more re- 
fined experimental setups. 
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